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Abstract 

It is shown that charged particle multiphcity distributions in restricted 
(pseudo)rapidity intervals in e~^e~ annihilation and in e~^p scattering at 
HERA are quite well described by the modified negative binomial distri- 
bution and its simple extension. 
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Recently it has been shown [1-5] that negatively charged particle multiplicity 
distributions in e+e~ annihilation are well described by the modified negative 
binomial distribution (MNBD) with the parameter A being practically energy 
independent and the parameter k approaching the value 7. It has been demon- 
strated also that the MNBD fairly well describes multiplicity distributions in 
lepton-nucleon interactions with a similar behaviour of the parameters [6]. The 
aim of this paper is to show that the MNBD and its simple extension describe 
quite well the charged particle multiplicity distributions in restricted rapidity 
intervals in e'^e" annihilation into hadrons [7-12], and charged particle multi- 
plicity distributions in restricted pseudorapidity intervals in e'^p scattering at 
HERA energies [13]. 

The MNBD is defined by the probability generating function 



where P„ is the probability to produce n particles. The MNBD can be con- 
sidered as a convolution of the Newton binomial (the numerator in (1)) and the 
negative binomial (the denominator in (1)). The probabilities Pn for the MNBD 
can be computed also using iteration relations given in [4] or the formulae given 
in [14]. We have checked that all three methods give the same P„ values within 
the computer precision limits. 

It has been assumed in [1,3,4] that the parameter k is the number of sources of 
particle production at some initial stage of the interaction; these sources develop 
independently of each other according to some branching process, producing in- 
termediate neutral clusters. In this model the parameter —A is equal to the 
cluster decay probability into a charged hadron pair (1 -|- A is the probability 
of cluster decay into pair of neutral hadrons). One can associate these sources 
with the widely discussed clans [15] or superclusters [16]. A priori it is not evi- 
dent that the same distribution should be valid for restricted rapidity intervals; 
indeed the initial sources can move with respect to each other and therefore 
contribute differently to the given rapidity interval, the same can happen for 
the intermediate neutral clusters. Let us assume that a fixed number of initial 
sources (less than or equal to k) contributes to a given rapidity interval (yi, 1/2), 
and both charged decay products of each neutral cluster hit the rapidity interval 
with the probability 62 and only one charged particle hits the interval with the 
probability ei. The probabilities ei and 62 are assumed to be the same for all 
neutral clusters, whatever their source. In this case the probability generating 




(1) 
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function for the charged particle muhiphcity distribution has the form (1), where 
X is replaced by 

¥P2(x) = l-ei(l-x)-e2(l-a;2). (2) 

The probability ei should increase with the size of the rapidity interval up to a 
maximum, and then go down to zero; whereas the probability €2 should increase 
monotonically from zero to one. One expects that this rapidity dependence can 
naturally explain odd-even fluctuations for rapidity intervals covering nearly 
the whole phase space. One expects also that the parameter k increases from 
some initial value (may be 1) to the maximum, observed over the whole phase 
space. For some rapidity intervals, when contributions from configurations with 
different number of sources are equally important, the above assumptions seem 
to be crude; in this case one can use a sum of distributions with different k 
values^. 

Analytical formulae for P„ can be obtained using the decomposition of the 
quadratic term in the the denominator of the probability generating function into 
the product of two linear terms. For k = 1 it leads to the following formulae for 
the probabilities: 

P„.l±^(!l±M, (3) 

1 + r (ei + ea) 



(l + r(ei + e2))(A-/32) 



I > 0, (4) 



where 




and 



pi = 'i+fi+h , P2 = 'j-\ri+h (5) 

^^ = l+r(e! + e2) ' = 1 + r (e' + e2) " 

For k greater than one we have calculated the P„ as a k-th convolution of the 
multiplicity distribution with k = 1, given by the formulae (3) and (4). It has 
been checked that for €2 = these calculations give the same results as the 
MNBD formulae. 

Probabilities P„ can be computed also using the recurrence relations for 
Gegenbauer polynomials C^{t) with the generating function 

00 

T.C::{t)z- = {l-tz + z')-' (7) 

n=0 



■^This approach, for example, is used in [17], where multiphcity distributions are 
parametrized by a weighted superposition of two negative binomials. 
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(see for example [18]). The case A = gives 



Pn+l = Pn + ^^±^6, P„_l (8) 

n + 1 n + 1 

with bi and 62 given by (6) and with initial conditions 

Po = (1 + r(ei + e2))-' , Pi = khP^ . (9) 

These recurrence relations can be used also for the case with r = if one 
replaces k by —k and r by A, and the final probabilities can be calculated 
by the convolution of these two distributions with A = and r = 0. We have 
checked that this method gives the same results as the method based on formula 
(4). 

The experimental data have been parametrized by the MNBD (e2 = 0) and 
its extension with €2 7^ 0. The parameter A was fixed at the value A = —0.76 
taken from fits of the e~^e~ data [4,5] and the parameter r was calculated from 
the mean charged multiplicity < n > using the relation 

r = A + , ^ , . (10) 

The integer parameter k has been tested in the interval from 1 to 9, and the 
parameters ei and €2 have been assumed to be nonnegative. One should note 
that the parameter A cannot be obtained directly from fits, indeed the same P„ 
can be obtained for other A values by simple rescaling of the ei, €2 and r. 

Tables 1 and 2 give the results of the fits to the charged particle multiplicity 
distributions in the symmetric rapidity intervals \y\ < Ucut for e^e~ annihilation 
for the TASSO [8], DELPHI [9] and ALEPH [11] dataQ. In some cases when the 
results with €2 7^ are not given, the fits have stopped on the boundary €2 = 0. 

One can see from tables 1 and 2 that the quality of the fits is quite qood 
for all data excluding the DELPHI measurement [9]. It is necessary to note 
that the x^/^DP values for the fits should be considered just indicative, since 
the full covariance matrix is not published and therefore the proper treatment 
of the correlations between measurements of the neighbour multiplicities is not 
possible]^. Nevertheless, for the TASSO measurements one can compare these 
values with the /NDF values obtained for negative binomial fits [8] made to 

The Pn for the TASSO measurements [8] were calculated from the KNO scaling function 
ip{z), compiled in the COMPAS data base, IHEP, Protvino. 

^Typical uncertainties (r.m.s.) in measured multiplicities in recent high energy experiments 
range from ~ 2 at rich = 8 to ^ 4 at rich = 30 [11]. 
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the unfolded distributions. The /NDF ratios for MNBD fits in general are 
smaller than for negative binomial fits. Especially it is true for the statistically 
most significant TASSO measurement at i/i = 34.8 GeV. The parameter k tends 
to rise with increasing size of the rapidity interval. The parameters ei and €2 in 
general follow the expected dependence on the size of the rapidity interval. 

The big /NDF ratios for the DELPHI measurements [9] are explained 
by the shoulder structure caused by the superposition of events with different 
number of jets (for the discussion see [17,19]). This shoulder structure, well 
described by the Monte Carlo program JETSET, is parametrized in [17] by a 
weighted superposition of two negative binomials. Following this approach we 
have parametrized multiplicity distributions in restricted rapidity intervals by 
the weighted superposition of two modified negative binomials. Let us denote 
the probabilities P„ for MNBD with the form (2), and £2 = by Pnik, A, ri, ei). 
Then the weighted sum is given by 

Pn = apnik, A, ri, ei) + (1 - a)pn{k + 1, A, r2, ei) , (11) 

where a is the probability for the contribution of the k sources with the pa- 
rameter ri, whereas 1 — a is the probability for the contribution of the k + 1 
sources with the parameter r2. The results of fitting the multiplicity distribu- 
tions measured by DELPHI [9] in symmetric rapidity intervals by the weighted 
sum (11) are given in table 3. One can see the 'x^ /NDF ratios are close to one; 
nevertheless it is necessary to note that the fits have four free parameters. We 
have tried also to fit by this parametrization the multiplicity distributions [9] in 
asymmetric rapidity intervals < y < Ucut ( or —ycut < Z/ < 0). The results are 
given m table 4. The x^/NDF ratios are near 1 for all intervals, excluding the 
intervals with t/cut = 1 and ycut = 1.5 : for these intervals, the ratios are near 2. 

The results of the fits of the MNBD, and its extension with 62 7^ 0, to charged 
particle multiplicity distributions in the single hemisphere for e~^e~ annihilation 
into hadrons [7-10,12] are given in table 5. The x^ /NDF ratios are of order 
1 for practically all measurements. One should note also that the fits to the 
DELPHI measurements [12] for hb events and all flavours have smaller x^ /NDF 
ratios than the "shifted" negative binomial fits given in [19]0. For example, the 
MNBD fit for the hh sample gives x^/NDF = 13.3/27, and the fit with €2 7^ 
gives x^ /NDF = 7.9/26, whereas the best negative binomial fit with a shift in 
multiphcity by 2 units has x^/NDF = 20/26 [19]. 

It is of interest to note that the k values for single hemispheres are close to 
the parameter k for MNBD fits to the multiplicity distributions in the whole 

^Thc shifts in multiplicity are due to the decays of the heavy 6-quarks. 
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phase space [4,5], whereas the parameter ei is rather close to the value 0.5. The 
second observation can be explained if one assumes that the intermediate neutral 
clusters hit one or other hemisphere with equal probability. 

The results of the fits to the HI data [13] in the pseudorapidity intervals 
1 < ?7 < r]cut are given in table 6. The x^/^DF values are in general close 
to 1 both for 62 = and for £2 7^ 0. It is of interest to note the the fits with 
€2 7^ have smaller /NDF than the MNBD fits mainly for the energy W 
above 150 GeV. 

In conclusion, it is shown that the MNBD and its simple extension describe 
multiplicity data in the restricted (pseudo)rapidity intervals for e+e~ annihi- 
lation and e^p deep inelastic scattering at HERA quite well. The DELPHI 
measurements [9] in symmetric and asymmetric rapidity intervals are fairly well 
parametrized by the weighted superposition of two modified negative binomial 
distributions with the parameters k and k+1. The dependence of the parameters 
on the size of the rapidity interval is explained by simple assumptions, extending 
the model proposed in [1,3,4]. Futher checks of this phenomenological model can 
be done, for example, by studying the multiplicity distributions for particles of 
the same charge. 
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Table 1: Results of the fits with fixed parameter A = —0.76 to the charged 
particle multiplicity distributions measured by TASSO [8] as a function of the 
centre of mass energy in restricted rapidity intervals \y\ < ijcut- The last 
column gives the /NDF for the negative binomial fits taken from [8]. 



yfs (GeV) 


Vcut 


k 


ei 


62 


XVNDF 


xVndf(nb) 




0.5 


3 


0.151±0.022 





12.4/14 


15/13 






2 


0.391±0.017 


0.106±0.029 


9.7/13 




14.0 


1.0 


5 


0.108±0.023 





37.2/20 


27/19 




1.5 


7 


0.164±0.023 





36.0/23 


33/22 






5 


0.372±0.045 


0.170±0.061 


36.8/21 






2.0 


9 


0.386±0.023 





85.0/25 


75/24 






5 


0.239±0.021 


0.544±0.018 


72.8/24 






0.5 


3 


0.181±0.022 





3.0/14 


4/13 






2 


0.411±0.020 


0.119 ±0.030 


4.8/13 




22.0 


1.0 


3 


0.421±0.027 





4.6/21 


19/20 






2 


0.470±0.045 


0.379±0.043 


11.7/20 






1.5 


5 


0.283±0.028 





6.5/26 


8/25 






5 


0.337±0.030 


0.054±0.91 


6.4/25 






2.0 


9 


0.186±0.025 





14.5/28 


13/27 






8 


0.269±0.057 


0.190±0.076 


11.9/27 






0.5 


3 


0.153±0.011 





10.1/15 


22/14 






2 


0.404±0.008 


0.109±0.015 


6.9/14 




34.8 


1.0 


3 


0.387±0.013 





49.3/26 


102/25 






2 


0.462±0.018 


0.403±0.017 


50.3/25 






1.5 


4 


0.415±0.015 





34.7/31 


74/30 






3 


0.406±0.034 


0.413±0.029 


30.7/29 






2.0 


8 


0.082±0.015 





28.1/34 


27/33 






6 


0.308±0.039 


0.204±0.048 


30.7/33 






0.5 


3 


0.080±0.022 





6.7/17 


8/16 






2 


0.377±0.015 


0.077±0.034 


7.1/16 




43.6 


1.0 


3 


0.315±0.026 





25.8/27 


44/26 






2 


0.451±0.035 


0.399±0.034 


21.4/26 






1.5 


5 


0.105±0.027 





45.6/32 


39/31 




2.0 


7 


0.117±0.027 





43.6/34 


39/33 
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Table 2: Results of the fits with fixed parameter A = —0.76 to the charged 
particle multiplicity distributions measured by DELPHI [9] and ALEPH [11] at 
^/s = 91.2 GeV in the symmetric rapidity intervals \y\ < Ucut- 



Collaboration 


Ucut 


u 

K 




^2 






0.5 


2 


0.155±0.018 





55.8/21 


DELPHI 


1.0 


2 


0.631±0.016 





111.3/35 




1.5 


3 
2 


0.481±0.014 
0.362±0.019 




0.730±0.018 


279.4/44 
295.6/43 




2.0 


4 
3 


0.425±0.014 

0.316±0.030 




0.606±0.027 


222.5/48 
218.9/17 




0.5 


2 


0.305±0.057 





5.1/19 


ALEPH 


1.0 


2 
2 


0.753±0.068 
0.552±0.145 




0.301±0.162 


11.6/33 
10.8/32 




1.5 


3 
2 


0.519±0.056 
0.418±0.089 




0.696±0.082 


32.7/39 
16.7/38 




2.0 


3 
2 


0.929±0.054 
0.335±0.066 




0.981±0.057 


10.5/45 
37.0/44 



9 



Table 3: Results of the fits by the superposition of the two weighted MNBDs 
to the charged particle multiplicity distributions measured by DELPHI [9] at 
= 91.2 GeV in restricted rapidity intervals \y\ < t/cut- 



Vcut 


k 


a 


ei 


n 


r2 


XVNDF 


0.5 


5 


0.735±0.022 


0.058±0.034 


6.87±4.62 


15.28±9.36 


11.0/18 


1.0 


5 


0.727±0.013 


0.140±0.024 


6.08±1.27 


12.31±2.29 


27.9/32 


1.5 


8 


0.681±0.010 


0.062±0.022 


6.50±1.24 


13.03±2.22 


51.1/41 


2.0 


9 


0.570±0.011 


0.131±0.024 


7.05±1.47 


13.19±24.92 


46.7/45 



Table 4: Results of the fits by the superposition of the two weighted MNBDs 
to the charged particle multiplicity distributions measured by DELPHI [9] at 
= 91.2 GeV in restricted rapidity intervals < y < i/cut- 



Vcut 


k 


a 


ei 


n 


r2 


xVndf 


0.5 


4 


0.790±0.025 


0.060±0.039 


3.96±2.24 


9.58±6.78 


3.6/12 


1.0 


4 


0.848±0.009 


0.061±0.029 


9.50±4.52 


23.83±11.38 


35.4/25 


1.5 


6 


0.769±0.008 


0.065±0.022 


8.05±2.98 


20.51±7.03 


55.8/29 


2.0 


9 


0.685±0.008 


0.055±0.021 


7.84±3.27 


19.51±7.48 


75.5/30 


3.0 


9 


0.423±0.017 


0.219±0.027 


2.24±0.41 


4.51 ±0.62 


36.5/30 


5.0 


9 


0.563±0.061 


0.467±0.025 


1.44±0.19 


2.01 ±0.16 


31.5/30 
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Table 5: Results of the fits with fixed parameter A = —0.76 to the charged 
particle multiplicity distributions in a single hemisphere for e'^e~ annihilation 
into hadrons. 



Experiment 


(GeV) 


k 


ei 


€2 


XVNDF 




14.0 


6 


0.566±0.009 





10.5/17 






4 


0.568±0.027 


0.292±0.024 


7.1/16 


TASSO [8] 


22.0 


8 


0.456±0.010 





5.5/17 






6 


0.535±0.024 


0.162±0.030 


5.0/16 




34.8 


8 


0.470±0.005 





26.4/23 






6 


0.541±0.015 


0.196±0.017 


24.7/22 




43.6 


9 


0.376±0.010 





10.6/23 






6 


0.483±0.033 


0.253±0.033 


12.5/22 


HRS [7] 


29.0 


9 


0.487±0.011 





17.8/18 






5 


0.576±0.047 


0.321±0.039 


24.5/17 


OPAL [10] with Ks, A 


91.2 


5 


0.653±0.018 





18.7/34 






5 


0.667±0.080 


0.245±0.078 


12.1/33 


OPAL [10] no Ks, A 


91.2 


5 


0.782±0.018 





7.6/32 






5 


0.684±0.086 


0.150 ±0.100 


7.0/31 


DELPHI [9] 


91.2 


6 


0.671±0.008 





98.7/34 






4 


0.309±0.023 


0.729±0.018 


80.7/33 






6 


0.844±0.011 





13.3/27 


DELPHI [12]. ^/^-ovoiits 


91.2 


5 


0.592±0.0(38 


0.152±0.062 


7.9/26 






6 


0.691±0.012 





24.9/29 


DELPHI [12], all flavours 


91.2 


5 


0.506±0.057 


0.412±0.048 


20.1/28 
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Table 6: Results of the fits with fixed parameter A = —0.76 to the charged par- 
ticle multiplicity distributions measured by HI [13] in restricted pseudorapidity 
intervals 1 < r} < rjcut- 



W(GeV) 


1] interval 


k 


ei 


£2 


XVNDF 


80^ 115 


1^2 


3 
1 


0.038±0.073 
0.547±0.015 




0.303±0.047 


0.8/13 
14.3/12 


1^3 


4 
3 


0.100±0.076 
0.339±0.102 




0.167 ±0.134 


1.8/17 
2.2/16 


1 H-4 


9 
5 


0.026±0.079 
0.317±0.215 




0.229±0.228 


0.5/18 
0.7/17 


1^5 


9 
3 


0.162±0.044 

0.346±0.056 




0.627±0.060 


1.6/18 

9.9/17 


115^ 150 


1^2 


2 
1 


0.199±0.035 
0.533±0.026 




0.285d=0.027 


6.9/14 
62.5/13 


1^3 


3 
1 


0.178±0.090 

0.588±0.060 




0.385±0.055 


2.8/18 

113.3/17 


1^4 


7 
7 


0.043±0.056 
0.174±0.048 




0.122±0.110 


4.5/20 
3.4/19 


lH-5 


9 
4 


0.154±0.047 
0.319±0.058 




0.515±0.056 


7.6/21 
12.8/20 


150^ 185 


1^2 


3 
3 


0.028±0.092 

0.182±0.055 




0.066±0.042 


2.9/14 

1.0/13 


1^3 


3 
2 


0.195±0.116 
0.442±0.170 




0.322±0.097 


4.0/20 
3.0/21 


1^4 


6 
6 


0.064±0.083 
0.162±0.083 




0.178±0.097 


4.4/22 
3.0/21 


1^5 


9 
9 


0.088±0.064 
0.218±0.034 




0.154±0.088 


3.4/22 
1.5/21 


185 ^ 220 


lH-2 


3 
3 


0.034±0.096 
0.173±0.064 




0.056±0.041 


2.4/14 
0.6/13 


1^3 


3 
2 


0.188±0.096 

0.366±0.075 




0.322±0.059 


3.8/21 

8.4/20 


1^4 


5 
4 


0.126±0.091 
0.232±0.060 




0.300d=0.067 


9.2/22 
8.1/21 


lH-5 


9 
5 


0.169±0.074 
0.233±0.060 




0.502±0.064 


8.7/23 
6.2/22 
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